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Abstract 

The density of orbital space debris constitutes an increasing environmental chal- 
lenge. There are three ways to alleviate the problem: debris mitigation, debris re- 
moval and collision avoidance. This paper addresses collision avoidance, by describ- 
ing a method that contributes to achieving a requisite increase in orbit prediction ac- 
curacy. Batch least-squares differential correction is applied to the publicly available 
two-line element (TLE) catalog of space objects. Using a high-precision numerical 
propagator, we fit an orbit to state vectors derived from successive TLEs. We then 
propagate the fitted orbit further forward in time. These predictions are compared 
to precision ephemeris data derived from the International Laser Ranging Service 
(ILRS) for several satellites, including objects in the congested sun-synchronous or- 
bital region. The method leads to a predicted range error that increases at a typical 
rate of 100 meters per day, approximately a 10-fold improvement over TLE's prop- 
agated with their associated analytic propagator (SGP4). Corresponding improve- 
ments for debris trajectories could potentially provide initial conjunction analysis 
sufficiently accurate for an operationally viable collision avoidance system. 

We discuss additional optimization and the computational requirements for ap- 
plying all-on-all conjunction analysis to the whole TLE catalog, present and near 
future. Finally, we outline a scheme for debris-debris collision avoidance that may 
become practicable given these developments. 
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1 Introduction 



Collisions in orbit pose a threat to spacecraft, to astronauts and to the 
global commons of near-earth space. Several collisions have already occurred 
between spacecraft and debris, while the Iridium 33/Cosmos 2251 collision 
of January 2009 represented the first documented satellite-satellite collision. 
Unfortunately, the cumulative number of collisions thus far is consistent with 



the prescient predictions of a runaway chain reaction (Kessler & Cour-Palais 



1978). 



Reducing debris from current and future space missions (debris mitiga- 
tion) is an important measure. However, even the unrealistic case of "no new 



launches" is inadequate to curb runaway growth (Liou & Johnson, 2008). 
There are two remaining ways to alleviate the problem: debris removal and 
collision avoidance. For the latter, one needs prediction precision to be suffi- 
ciently accurate so as not to imply an unwieldy number of collision avoidance 
maneuvers. In addition, one needs a means of doing debris-debris collision 
avoidance since applying collision avoidance only to the subset of conjunc- 
tions involving maneuverable spacecraft does not suffice to curb the growth. 
This paper describes a method to increase the orbit prediction accuracy based 
on publicly available TLEs. 



Many satellite owner operations have inadequate (if any) processes for con- 
junction assessment and collision avoidance, since they would have to screen 
their asset (s) against all other space objects. The only source of knowledge 
at their disposal for the majority of other objects is the publicly available 
two-line element (TLE) sets. But predictions based on TLEs using the as- 
sociated analytic propagator (SGP4) are not sufficiently accurate to warrant 
maneuvering to avoid potential collisions: they imply an unacceptably large 
number of potential collisions per space object, each of which has very low 
probability. The problem is similar for debris-debris conjunctions except then 
both objects, not just one, are subject to these imprecisions. 

To address this TLE/SGP4 accuracy problem, we investigated several meth- 
ods to improve the propagation errors for non-maneuvering orbital objects 
whilst using only TLEs as input data. The following research target was set: 
"increase the predictive accuracy for orbital objects, using only historical TLE 
data, such that it enables operational conjunction assessment for collision 
avoidance" . 



In the following sections, we describe one approach to this target, assess its 
accuracy, and discuss the requirements for extension to the entire space object 
catalog. Finally, we propose a new method of debris-debris collision avoidance 
enabled by long-term high-accuracy conjunction assessment. 
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2 Method: TLE orbit fitting and propagation 



ESA uses TLEs from the publicly available catalog to initially screen their 
sun-synchronous orbit (SSO) spacecraft ERS2 and ENVISAT for conjunctions 
( Flohrer, Krag, &: Klinkrad| 2009). Telemetry from their operational space- 
craft provide precision orbital ephemeredes (POEs) for those spacecraft. For 
screening against all other potentially conjuncting objects, only TLEs are 



used. Flohrer, Krag, & Klinkrad (2008) describe a method to estimate er- 



ror covariances of TLEs in order to quantify collision probability assessments. 
Their method provided inspiration for the present work. Here we extend and 
adapt their approach: based solely on the object's historical TLEs, we improve 
the accuracy of the object's predicted position, as opposed to quantifying the 
accuracy of the object's SGP4 propagation errors. 

Our method is essentially to use TLE data as "pseudo-observations" and 
to fit an orbit to these pseudo-observations using a high-precision special per- 
turbations propagator and traditional batch least-squares differential correc- 
tion. The fitted orbit is then propagated into the future using the same high- 
precision orbit propagator. The prediction accuracy is assessed by comparison 
with precision orbital ephemeris (POE) data from the International Laser 



Ranging Service (Pearlman, Degnan, & Bosworth, 2002). 



For each object we wish to analyze we choose a time window with two 
sections: a fitting period and a subsequent prediction period. The fitting pe- 
riod is initially set to ten days, typical of the period over which U.S. Space 
Surveillance Network observations are fitted when generating TLEs for LEO 
objects (Danielson et al. 2000). Section 4 details a more principled approach 



to determining the fitting period. The prediction period is 30 days. 

For a given object, all TLEs with epochs within its entire window are ob- 
tained. However, only TLEs with epochs within the fitting period are used for 
the fitting process; those in the prediction period are only used for validation 
and test. 



Using all TLEs from the fitting period and interpolating using SGP4, we 
generate a series of order 100 state vectors ("pseudo-observations") equally 
spaced in time within the fitting period. We initialize the differential cor- 
rection with a state vector derived from the first TLE in the fitting period, 
and then propagate that using a high precision propagator. We apply batch 
least-squares differential correction to minimize the RMS error in the radial, 
in track and cross track (RIC) components of relative positions between the 
fitted orbit and the pseudo-observations. The trajectory obtained from the 
converged differential corrector is our fit. 

The high precision propagator incorporates a 60 x 60 EGM2008 gravity 
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field including solid Earth tides, point masses for the solar and lunar gravity 
fields and a NRL-MSISE 1990 atmospheric drag model. We assume constant 
values for solar F10.7 and Ap equal to their averages over the fitting period. 
Solar radiation pressure is modelled using a simple the biconic approximation 
with Earth and Moon as eclipsing bodies. 



3 Results: Prediction Accuracy 



We applied our method to four non-manoeuvring spacecraft for which POE 



data were readily available (Ries, 2009): Stella, Starlette, Ajisai and Etalon-2. 



The summary of basic orbital and physical properties for these spacecraft are 
shown in Table 1. They were picked in part for their range of orbital altitudes. 

Table 1 

Key Parameters of Geodesic Spacecraft Used 
Satellite Perigee Eccentricity Period Inclination Mass Diameter 
(km) (min) (deg) (kg) (m) 



Stella 


800 


0.0206 


98.6 


101 


48 


0.24 


Starlette 


812 


0.0206 


104 


49.8 


47 


0.24 


Ajisai 


1490 


0.0010 


116 


50.0 


685 


2.15 


Etalon 2 


19120 


0.0007 


675 


65.5 


1415 


1.29 



Figure 1 shows examples for two of these spacecraft - the highest and lowest 
altitude - comparing our method's predictions to those of SGP4. The TLEs for 
these objects exhibit quite different behaviours. Nevertheless, in both cases the 
TLEs propagated with SGP4 (black) depart rapidly from the true position of 
the satellite whereas our numerically fitted and numerically propagated orbits 
(red) maintain higher accuracy, particularly over the long term. 

Precision orbit ephemeredes ("truth data") define the x-axis for plotting, 
but truth data are not used in the fitting process. Nor are they used for tuning 
hyperparameters (e.g. length of fitting interval, see Section 4). Improvements 
in prediction accuracy are typically seen in all three axis (RIC), especially 
more than day into the future. For predictions of less than one day the errors 
appear to be dominated by TLE bias (with respect to truth) which our method 
cannot remove. 

For each of the four satellites, we performed 50 runs of our method using 
different starting dates distributed uniformly in 2004. Summary statistics of 
these runs using "box- and- whisker" plots appear in Figures 2 and 3. The 
TLE+SGP4 errors grow more rapidly in every case. The typical prediction 
errors for our method are 3km at 30 days out, corresponding to ~ 100 m/day 
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Fig. 1. Typical position errors with respect to POE ("truth") using our method 
(red), and using SGP4 (black) for Stella (upper panels) and Etalon-2 (lower panels). 
The prediction period is < t < 30 days. The fitting period is — 10 < t < days 
for Stella and —90 < t < days for Etalon-2. All fitting uses only publicly available 
TLEs. Truth data are only used for plotting. Selection of the length of the fitting 
period is discussed in Section 4. Also shown are the updated TLEs (no prediction) 
in blue. 

prediction error growth. The worst case is Stella where the median error grows 
< 150 m/day and the best is Etalon at < 30 m/day. The maximum error for 
our method over all four satellites is < 300 m/day. These compare to typical 
growth of 1,500 m/day for TLE's propagated with SGP4 for these objects. 
The ratio of median prediction error using our method vs. SGP4 at 30 days 
range from 3 for Etalon-2 to 50 for Starlette, and averages 15. 

Thus for these satellites, our method exhibits approximately one order of 
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Fig. 2. "Box and whisker" plots of range errors for predictions using our fit method 
(red), and using SGP4 (black) for four satellites. Each panel summarizes 50 runs 
corresponding to 50 random starting dates in 2004. The x-axis is prediction time in 
days. The y-axis is prediction error in km. The lower and upper whiskers extending 
from each box bound the minimum and maximum prediction errors respectively. 
The lower and upper box edges bound the 25 th and 75 th percentiles, respectively. 
The glyphs inside the boxes mark the median prediction error. All errors are plotted 
with respect to "truth" (i.e. POEs obtained from the ILRS). Truth data are only 
used for plotting. All fitting uses only publicly available TLEs. 

magnitude improvement in prediction accuracy over TLEs propagated with 
SGP4. Since the improvement is in all three directions (RIC) the resultant 
decrease in position covariance ellipsoid volume is likely to result in 2 to 3 
orders of magnitude reduction in false positive conjunctions. 

Since the TLEs for these objects have an instantaneous range bias from 
truth of 0.8 ± 0.3 km the error growth of approximately 1.5 km/day means 
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Fig. 3. Additional summaries of range errors for predictions using our method (red), 
and using SGP4 (black) for four satellites. On the LHS, just our method is plotted 
for a full 30 days and on the RHS both our method and the SGP4 propagated TLEs 
are shown but just for the first 10 days (5 in the case of Starlette since the TLE was 
diverging rapidly) . Each box and whisker plot summarizes 50 runs corresponding to 
50 random starting dates in 1994. The x-axis is prediction time in days. The y-axis 
is prediction error in km. All errors are plotted with respect to "truth" (i.e. POE 
obtained from the ILRS). Truth data are only used for plotting. All fitting uses only 
publicly available TLEs. 



that only after approximately one day can one detect the benefits of our 
method. Our fitting method has a similar initial range bias to TLEs but error 
growth of only 100 m/day. These errors are similar to those published in the 
open literature evaluating the high accuracy special perturbations catalog(s) 



7 



maintained by the US Space Command (Neal, Coffey, & Knowfes, 1997 Coffey 



et al. 1998 Boers et al. 2000). Table 2 compares the accuracy of various 



tracking data and prediction methods. 
Table 2 



Accuracies of Different Sensors and Prediction 


Methodologies. 




Method 


Prediction 


Accuracy 


Accuracy 




(days) 


(m) 


(m/day) 


Laser Ranging (ILRS "truth" ) 1 





0.1 


N.A. 


Fence (raw direction cosines) 2 





10 


N.A. 


High Accuracy catalog + SP 3 


10 


500-2000 


50-200 


TLEs + SGP4 4 


10 


1,000-30,000 


100-3,000 


TLEs + new scheme 


10 


500-2,000 


50-200 



4 Fitting Debris objects 



The majority of tracked space objects (i.e debris) do not have precision 
orbital ephemerides - they are not tracked by the ILRS and they transmit 
nothing. Many are, however, tracked by the United States Space Surveillance 
Network (SSN). Unfortunately, the only publicly available source of data de- 
rived from these SSN observations are, at present, the TLEs accessible at 
space-track.org (and archived e.g. at celestrak.com). 

While there is currently no publicly available source of high-accuracy tra- 
jectory for debris objects, we believe the method described above could be 
applied to refine the predictions of debris objects' future positions based only 
upon comparatively inaccurate TLEs, facilitating conjunction assessment and 
collision avoidance. However, there are several issues which complicate the 
application of our procedure to debris: 

The first set of issues arise when TLEs are the only source of data - we have 
no "truth" in order to perform validation and test. This is solved by the use 
of TLEs from beyond the fitting period as a proxy for truth data. 



Our initial manual experiments suggested that different fitting periods were 
optimal for different orbit categories. These results are consistent with Alfriend 
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Fig. 4. TLEs from beyond the fitting period can be used instead of precision 
ephemeris (POE) data when optimizing the length of the fitting period. Left: predic- 
tion error with respect to future TLEs ('"proxy") vs. fitting period. Right: prediction 
error with respect to POE data ("truth") vs. fitting period. Each box and whisker 
plot summarises 50 runs corresponding to 50 random starting dates in 1994. 



et al. (2002); Danielson et al. (2000). We made a more principled investigation, 
utilizing cross-validation, to optimize the length of the fitting interval on a per- 
object basis using only TLEs for both fitting and validation. Figure 4 shows 
how a fitting period of 10 days was determined to be optimal for satellite 
Ajisai. The important point shown in the figure is that a 10 day fitting period is 
optimal both when validating predictions with truth data and when validating 
predictions with future TLEs. And the 10 day fitting period is optimal whether 
predicting 5, 10 or 30 days into the future. Thus it appears "future" TLEs can 
be used as a proxy for truth data when optimizing the length of the fitting 
period on a per object basis. 

The second set of issues is that for debris objects one must also solve for 
area to mass ratio (drag) and validate this as well using only future TLEs. We 
are currently performing these experiments. 
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5 Outlook 



5. 1 Methodological Improvements 



Further accuracy improvements could come from a principled analysis of 
the number of fitting points and their weighting as a function of location 
within the fitting period, the removal of outlier data (i.e. filtering) and other 
techniques from statistical orbit determination, treating TLEs as observations. 



Also, recent work Legendre et al. (2008); Muldoon & Elkaim (2008) suggests 



that there may be simple global transformations of TLE data that can improve 
TLE/SGP4 accuracy. These systematic corrections might also be further re- 
fined on the basis of orbital parameters, resulting in a look-up-table used to 
correct TLEs as a function of their orbital regimes. 



In addition, the scheme could be tested on actual known conjunctions in 
order to analyse its efficiency for conjunction assessment. These issues will be 
the subject of a future paper. 



5.2 Computational Scaling 



As an initial assessment of the computational requirements needed to per- 
form conjunction analysis for all space objects in the publicly available catalog, 
a simple conjunction analysis system was parallelized/ 1 "] on the Pleiades super- 
computer at NASA Ames Research Center. Using approximately 200 CPU 
cores, checking all objects against all objects for conjunctions 7 days into the 
future takes about a minute. The conjunctions involving spacecraft precisely 
match the results from SOCRATES (Kelso & Alfano, 2005). The same system 
takes about 40 minutes to perform 7 days of all-on-all conjunction assessment 
of a "simulated S-band fence" catalog (obtained from the NASA Orbital De- 
bris Office) containing "pseudo-TLEs" for approximately 2.5 million objects 
of size > 2cm. This demonstrates the feasibility of doing all-on-all conjunc- 
tion assessment on the present and expected near-future catalogs. Scaling the 
fitting method discussed above to the future catalog will be the subject of a 
later paper. 



1 This work was performed by Chris Henze of the NASA Advanced Supercom- 
puting Division at NASA Ames Research Center. Details will be forthcoming in a 
future publication. 
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5.3 A New Scheme for Debris-Debris collision avoidance 



An along track AV of 1 cm/s provides of order 1km displacement per day. 
Since error growth for our method is 100m per day, such small manoeuvres 
might suffice for collision avoidance. For manoeuvres this small one could use 
radiation pressure to impart the necessary momentum. To verify the conjec- 
ture, we increased the flux of radiation by a factor of lOx the solar radiation 
constant for 10 minutes on a typical debris object (area to mass ratio 0.1 
m 2 /kg) in SSO. We then propagated that for a further 5 days. The displace- 
ment was > 0.1 km/day along track compared to a reference object that was 
not so irradiated. Two cumulative such events, if one ensured appropriate 
geometry, would thus lead to the necessary along-track displacement p~| 

A factor of lOx the solar rational constant could conceivably be applied from 
the ground. For example, a lOkW class laser attached to a lm class telescope 
with appropriate tracking, could, (considering diffraction limits only) provide 
the necessary radiation pressure for a debris object of area <lm 2 in SSO. One 
would need to perform a detailed engineering analysis of the feasibility of such 
a device, particularly the viability of maintaining small beam divergence and 
tracking. 

Debris-debris collision avoidance could eliminate the need for an active de- 



bris removal program. The results of Liou & Johnson (2009) imply, to first 



order, that performing debris-debris collision avoidance for 5 additional care- 
fully selected objects per annum would curb debris field growth. Furthermore, 
any additional debris-debris collision avoidance would start reducing the net 
debris density, eventually to below a critical level, because debris creating col- 
lisions are being nulled whereas atmospheric drag continues to cause objects 
re-entry. Thus this method could provide not just a stop gap measure, but a 
permanent solution. 



6 Conclusions 



This fitting and propagation method, based solely on TLEs, could poten- 
tially provide initial conjunction analysis sufficiently accurate for an opera- 
tionally viable collision avoidance system. If similar improvements are possible 
for debris objects, as we expect, then it could be used for all-on-all conjunction 
assessment. However, to ensure that the scheme works in an operational set- 
ting, it would be necessary to have additional data: the so-called "uncorrelated 

2 This idea was first conceived in discussions with Dr. Riidiger Jehn from ESA, 
July 2009. 
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objects" (a.k.a. "analysts set') which account for approximately 30 percent of 



potential conjunctions in LEO (Neuman, 2008). Further, we show that all-on 



all conjunction assessment is possible with moderate computer infrastructure, 
even with the large increase in size of the catalog of tracked objects that is 
expected in the next few years. 



Finally, given accurate predictions, we claim debris- debris collision avoid- 
ance may be possible by externally inducing small manoeuvres using radiation 
pressure from 10 kW class power density ground-based lasers. If feasible, this 
could negate the need for a large scale debris removal program. 
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